DESY 01-130 
TTP01-21 
|hep-ph/0109084 
September 2001 



43 
Oh 



X 



Determination of a s and heavy quark masses from 
recent measurements of R(s 



J.H. Kiihn° and M. Steinhauser 6 

O 

(N 

£ 

o 



(a) Institut fiir Theoretische Teilchenphysik, 
Universitat Karlsruhe, D- 76128 Karlsruhe, Germany 



CT) ■ (b) II. Institut fiir Theoretische Physik, 

^-j- | Universitat Hamburg, D-22761 Hamburg, Germany 

00 ' 
O 

ON . 

O Abstract 



In this paper we compare recent experimental data for the total cross section 

a(e + e~ — * hadrons) with the up-to-date theoretical prediction of perturbative QCD 
for those energies where perturbation theory is reliable. The excellent agreement 
suggests the determination of the strong coupling a s from the measurements in 
43 ■ the continuum. The precise data from the charm threshold region, when combined 

with the recent evaluation of moments with three loop accurracy, lead to a direct 
determination of the short distance MS charm quark mass. Our result for the strong 
coupling constant af ] (5 GeV) = 0.235±g;gJ? corresponds to af ] {M z ) = 0.124t^, 
for the charmed quark mass we find m c (m c ) = 1.304(27). Applying the same ap- 
proach to the bottom quark we obtain mb(rrih) = 4.191(51) GeV. Whereas our result 
for a s (Mz) serves as a useful cross check for other more precise determinations, our 
values for the charm and bottom quark masses are more accurate than other recent 
analyses. 

PACS numbers: 12.38.-t 14.65.Dw 14.65.Fy 



1 Introduction 



The strong coupling constant and the quark masses are the basic input parameters for 
the theory of strong interaction. Quark masses are an essential input for the evaluation 
of weak decay rates of heavy mesons and for quarkonium spectroscopy. The prediction of 
these masses is an important task for all variants of Grand Unified Theories. To obtain 
the values in a consistent way from different experimental investigations is thus a must 
for current phenomenology. 

During the past years new and more precise data for a(e + e~ — > hadrons) have become 
available in the low energy region between 2 and 10 GeV. At the same time increasingly 
precise calculations have been performed in the framework of perturbative QCD (pQCD), 
both for the cross section as a function of the center-of-mass energy y/s, including quark 
mass effects, and for its moments which allow for a precise determination of the quark 
mass. A fresh look at the determination of the strong coupling and the quark masses 
based on these new developments is thus appropriate. 

For the determination of a s we shall use energy regions sufficiently far away from the 
threshold where pQCD plus local duality is expected to give reliable predictions. For the 
determination of the quark masses we will use sum rules where the input depends heavily 
on the threshold region with its rapidly varying cross section. 

Sum rules as tool to determine the charm quark mass have been suggested by the 
ITEP group long ago |[|. Subsequently these methods have been developed further 
and frequently applied to the bottom quark. Most of these later analyses concentrated 
on using relatively high moments which are less sensitive to the continuum contribution 
and exhibit a very strong quark mass dependence. However, this approach requires the 
proper treatment of the threshold, in part the resummation of the higher order terms of 
the Coulombic binding and a definition of the quark mass adopted to this situation like 
the potential- or lS'-mass |3|, |4j. Recently improved experimental results in the charm 
threshold region have become available ||, and at the same time the moments have 
been evaluated up to order a 2 s 0]. A fresh look at the evaluation of the charm quark 
mass with the help of sum rules is thus an obvious task. We will concentrate on low 
moments, say up to the fourth one and present the results for the fifth to the eighth 
moment just for illustration. This is a natural route to determine directly a short distance 
mass, say m c (m c ) or even m c (3 GeV) in the MS scheme as advocated in the original 
papers ^ In the present work we concentrate on the mass in the MS scheme, m c (/i), 
and adopt fi = 3 GeV as our default value, a scale characteristic for the present problem 
and sufficiently high to ensure convergence of the perturbative series. As already discussed 
in 0, fixed order pQCD is adequate, if only low moments are used in the analysis. 

Recently the determination of the charm quark mass got quite some attention. In H] 
the pole mass M c has been determined from the comparison of the direct determination 
of the hadronic contribution to Aa(Mz) with the determination using analytical contin- 
uation. This leads to the range M c = 1.33 — 1.40 GeV 0. The pole mass has also been 
determined in [TO] using QCD sum rules in combination with nonrelativistic QCD. Their 
result reads M c = 1.70(13) GeV and is significantly higher than the one of Ref. ||. The 
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MS mass given in [[HJ reads m c (m c ) = 1.23(9) GeV. On the basis of so-called Chauchy 
sum rules, recent experimental data and analytical results for the three-loop photon po- 
larization function the MS charm quark mass has been determined in JTI] with the result 
m c {m c ) = 1.37(9) GeV. In |L2| pseudo-scalar sum rules have been used to determine si- 



multaneously the decay constant fc and the charm quark mass. For the latter the value 
Tn c (m c ) = 1.10(4) GeV is given in the abstract of [JTJ] which is significantly lower than 
the other evaluations. 

We also want to mention that there is a recent lattice evaluation of the charm quark 
mass [13[] with the result m c (m c ) = 1.26(4) (12) GeV. The first error corresponds to the 
statistical and the second one to the systematical uncertainty. Note that this result is 
derived in quenched QCD and the corresponding uncertainty is not included. 

In the remaining part of the Introduction we would like to fix the notation and define 
the quantities we are dealing with in the remainder of the paper. 

It is convenient to normalize the radiatively corrected hadronic cross section and to 
define the ratio 

„. , a(e + e~ — > hadrons) . . 

R(s) = (1) 

Opt 

where a pt = 47ra 2 /(3s). As an inclusive quantity R{s) is conveniently obtained via the 
optical theorem from the imaginary part of the polarization function of two vector currents 
via 

R(s) = 127rlm[n(g 2 = s + ie)] , (2) 
where IT(g 2 ) is defined through 

-9V + ») n(g 2 ) = iJdxe"*(0\Tj r (x)ji(0)\0), (3) 



( 



with being the electromagnetic current. 

The perturbative expansion of R(s) can be written as 



fl(.) = E < + ^< + (^) 2 fl g) + f^l) S < + (4) 



7T V 7T / V 7T 



where the summation is performed over all active quark flavours Q. (We ignore the small 
singlet contribution at order a^.) For a comprehensive compilation of the individual pieces 
we refer to jnj |15| where also explicit results are given. We want to stress that the full 
quark mass dependence is available up to order a 2 s @, . In the case of Rq the first three 
terms in the high-energy expansion are known [TO . 

Our theoretical predictions are based on Eq. @ where the up, down and strange quark 
masses are taken to be massless and for the charm and bottom quark the respective pole 
masses are chosen as input. If not stated otherwise we will use the following input values 
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for the evaluation of R(s) 



af\M z ) = 0.118 ±0.003, 

M c = (1.65 ±0.15) GeV, 

M b = (4.75 ±0.20) GeV, (5) 

which cover the full range of all currently accepted results. 

At several places of our analysis the renormalization group functions and the matching 
conditions for the masses and the strong coupling are needed in order to get relations 
between different energy scales. The corresponding calculations are performed using the 
package RunDec [ |I~7j . 



The outline of the paper is as follows: in Section |2] we compare the experimental data 
of R(s) with the theoretical prediction and determine a s from the continuum data below 
y/s = 3.73 GeV and above y/s = 4.8 GeV. The measurement in the charm threshold region 
is used for the determination of the charm quark mass in Section [5| Similar considerations 
are used in Section [| in order to obtain the bottom quark mass. Section || contains our 
conclusions. 



2 The continuum region 

As stated in the Introduction, we distinguish two energy regions: first, the continuum 
region where pQCD and local duality are expected to give reliable predictions for the 
hadronic cross section and, second, the charm threshold region starting from the D meson 
threshold at 3.73 GeV up to approximately 5 GeV, where the cross section exhibits rapid 
variations, plus the J/\V and \V' resonances. The former will be mainly sensitive to the 
value of a s . The latter will be used to evaluate moments and to determine the charm 
quark mass. For the present analysis the continuum region covers the BES data points 
from 2 GeV up to 3.73 GeV and the data from BES §, MD-1 and CLEO |TJ between 
4.8 GeV and 10.52 GeVQ. As is evident from Fig. | pQCD with af\M z ) = 0.118 provides 
an excellent description of all recent results. 

Below 3.73 GeV only u, d and s quarks are produced and the 0(a z s ) approximation 
for massless quarks is adequate for a description of R. The effective number of flavours 
is chosen to be = 3 and virtual charm quark effects are taken into account (for a 
compilation of the relevant formulae see Ref. |jT4], |T5|). Above charm threshold u, d and 



s quark production is calculated as before, however, with n/ = 4. Up to 0{a 2 s ) the 
prediction for the charm quark production incorporates the full M c dependence. Starting 
from order a 2 s also the M c dependence of "secondary" charm production has to be taken 
into account. This includes diagrams of the type in Fig. ^|(a) as well as those from 
Fig. ^](b). In addition we include 0(a^) terms from the expansion in (M 2 /s) n with n = 0, 1 

1 We limit this analysis to more recent results from BES [pj, those from MD-1 |l8j and from CLEO |l9| 
with systematic errors of typically 4.3%, 4% and 2%, respectively. Older measurements, in particular 
those from SPEAR and DORIS, are consistent with the new results. However, with their significantly 
larger errors they do not provide additional information. 
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Figure 1: R{s) in the energy range between 1.8 GeV and 11.0 GeV. The solid line corre- 
sponds to the theoretical predictions adopting our central values for the input parameters. 
The theoretical uncertainties are indicated by the dashed curves which are obtained from 
the variation of the input parameters as described in the text. The two error bars on the 
data points indicate the statistical (inner) and systematical (outer) uncertainty. 

u,d,s 





Figure 2: Feynman diagrams contributing to R(s) at order a 2 . A secondary charm quark 
pair is produced through gluon splitting. 



and 2 [14], 16|. Last not least contributions from virtual c quarks {^/s < 3.73 GeV) and b 
quarks (-^i < 10.52 GeV) have been calculated which are suppressed ~ (a s /7r) 2 s/(4M 2 ) 
and decouple for s AM 2 . These are included in the fifth column of Tab. |l[ Pure 
QED final state radiation is tiny and taken into account for completeness. (For a related 
analysis at yfs = 10.5 GeV see |2(J.) 

The magnitude of the most important terms is shown for several characteristic ener- 
gies in Tab. [1] and compared to the experimental results. All numbers are given for the 
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(GeV) 


n f 




u, d, s 


c, b 


R th (s) 


R exp (s) 


2.000 


3 


0.297 


2.202 


0.002 


2.203 


2.180± 0.070 


3.730 


3 


0.226 


2.154 


0.002 


2.156 


2.100± 0.080 


4.800 


4 


0.215 


2.146 


1.607 


3.753 


3.660± 0.140 


8.900 


4 


0.180 


2.122 


1.440 


3.562 


3.579± 0.066 


10.520 


4 


0.172 


2.117 


1.427 


3.544 


3.560± 0.010 



Table 1: In the last two columns the experimental values for R(s) are compared with the 
theoretical prediction, R th (s), for a selected number of energies where fi 2 = s has been 
adopted (see text). For completeness we list the number of active flavours, rif, and the 

values of aS 1 (y/s). 



reference values as specified in Eq. (Q). The band in Fig. [l| is obtained from the corre- 
sponding errors and the theoretical uncertainty from the variation of the renormalization 
scale between fi 2 = s/4 and /i 2 = 4s which have been added quadratically. The excellent 
agreement between prediction and measurement suggests that one might determine a s 
from this analysis. 

Assigning a fully correlated systematic error of 4.3% to the BES data |21| we find 

«i 3 >(3GeV) = 0.3691^;^, (6) 

from the BES data below 3.73 GeV and 

a«(4.8GeV) = 0.183lSt°;°f r , (7) 

from the point at 4.8 GeV. The uncertainties refer to the uncorrelated (statistical) and 
correlated errors, respectively. The same analysis is applicable to the MD-1 results. On 
the basis of the R values listed for six energies in their table 4 and assuming fully correlated 
systematic errors we find 



a 



( 4 )(8.9GeV) = 0.201±g^i±8:a 



From their average value i?((8.9 GeV) 2 ) = 3.578±0.021±0.140 as quoted in their summary 
we obtain^ 



a 



( 4 ) i 



8.9 GeV) = 0.193j2$S82£, (9) 



consistent with Eq. (§) but with smaller errors. From the CLEO value R( (10.52 GeV) 2 ) = 
3.56 ± 0.01 ± 0.07 we deduce 

a?) (10.52 GeV) = 0.186^^ ■ (10) 



2 This result disagrees from their value a{ 4 \8.9 GeV) = 0.174 ± 0.039 as a consequence of their 
inadequate treatment of the quark mass effects. 
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We combine the results of Eqs. (|7|), (|9]) and flTOD by assuming uncorrelated sys- 
tematic errors and evolving them to a common scale of 5 GeV. Combining the errors in 
quadrature we find 



c4 4) (5 GeV) = 0.235^ 



0.047 • 



Evolving this value up to Mz 



af\M z 



n 1 94.+ - 011 

u - iZ4 -0.014 > 



121 



good agreement with other determinations 22, 23 is observed. 



3 Charm quark mass determination from the thresh- 
old region 

The first method we want to use for the determination of the charm quark mass is based 
on the direct comparison of theoretical and experimental moments of the charm quark 
contribution to the photon polarization function as defined in Eq. (Q). In the limit of 
small momentum the latter can be cast into the form pj 

n c (g 2 ) = Q 2 c ^- 2 Y,CnZ n , (13) 
167r ^0 

with Q c = 2/3 and z = q 2 /{^m 2 ) where m c = m c (/x) is the MS charm quark mass 
at the scale \x. The perturbative series for the coefficients C n up to n = 8 is known 
analytically ||, [jj up to order a 2 . The coefficients also depend on the charm quark mass 
through logarithms of the form l mc = ln(m^(/i)//i 2 ) and can be written as 



(14) 



In Tab. |2| the individual coefficients are given in numerical form. They essentially consti- 
tute our theoretical input. We define the moments 



r 2 / i \ n 



(15) 

q 2 =0 



which leads to 
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71 


i 
1 


9 


q 
o 


A 


E 




D 


7 
1 


Q 
O 


MO) 


1.0667 


0.4571 


0.2709 


0.1847 


0.1364 


0.1061 


0.0856 


0.0709 


r-(io) 


2.5547 


1.1096 


0.5194 


0.2031 


0.0106 


-0.1158 


-0.2033 


-0.2660 


M 11 ) 


2.1333 


1.8286 


1.6254 


1.4776 


1.3640 


1.2730 


1.1982 


1.1351 


M^o) 


2.4967 


2.7770 


1.6388 


0.7956 


0.2781 


0.0070 


-0.0860 


-0.0496 


rv(2i) 


3.3130 


5.1489 


4.7207 


3.6440 


2.3385 


0.9553 


-0.4423 


-1.8261 


^V(22) 


-0.0889 


1.7524 


3.1831 


4.3713 


5.3990 


6.3121 


7.1390 


7.8984 



Table 2: Coefficients of the photon polarization function in the MS scheme as defined in 
Eqs. (|13|) and fll4|). n/ = 4 has been adopted which is appropriate for the charm threshold. 



With the help of a dispersion relation we establish the connection between the polar- 
ization function and the experimentally accessible cross section R c (s). In the MS scheme 

n ^ 2 ) = wif^-r^ + ^ih ^ ( 17 ) 

127T J s(s — q z ) lD7r^ 
which allows to determine the experimental moments 

= J^lMs)- (18) 

Note, that the last term in Eq. fllTD which defines the renormalization scheme disappears 
after taking derivatives with respect to q 2 . Equating Eqs. (|l"6|) and (18) leads to an 
expression from which the charm quark mass can be obtained: 

1 ( C \ 1 ^ 2n " > 
mM = - 2 {j^) ■ (19) 

As a second method we consider the ratio of two successive moments which leads to 



m c (/i) -- 



Mf p C n+l (2Q) 



Here the normalization uncertainty of the experimental data is largely cancelled. Both in 
Eqs. fll"9f) and fl20"D one has to remember the m c dependence of C n . 

We have checked that the nonperturbative contribution from gluon condensates [I], 
can be neglected within the present accuracy. 

Let us now turn to the extraction of the experimental values for the three different 
contributions which enter the right-hand side of Eq. (0): the J/\l/ and resonances 
(M C n P,ms ), the charm threshold region between 2M Do ps 3.73 GeV and ^fs[ = 4.8 GeV as 
measured by the BES experiment [Q (.M° xp ' cc ), and the continuum contribution above si 
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The resonances are treated in the narrow width approximation which corresponds to 
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JP-(.) = ° J S(s-Ml), (21) 

cr \a(s) 



with a 1 = 137.0359895. For both resonances we use (a/a(s)) 2 ~ 0.9562 and 
Mj/* = 3.09687(4) MeV, F J J^ = 5.26(37) keV, M 9 > = 3.68596(9) MeV and F? = 
2.12(18) keV 

In the charm threshold region we have to identify the contribution from the charm 
quark, i.e. we have to subtract the parts arising from the light u, d and s quark from the 
data. Technically this is done by determining a mean value for R u d s = Ru + Rd + Rs from 
the comparison of theoretical predictions and the BES data between 2 GeV and 3.73 GeV 
and using the theoretically predicted energy dependence to extrapolate into the region 
between 3.73 GeV and 4.8 GeV ||25|| . This value is subtracted from the data before the 
integration is performed. Alternatively, one could adopt the massless prediction for R(s) 
up to order without taking into account the data below y/s = 3.73 GeV. We checked 
that both approaches lead to the same final result. 

In the continuum region above y/s = 4.8 GeV there is only sparse and quite unprecise 
data. On the other hand pQCD provides reliable predictions for R(s), which is essentially 
due to the knowledge of the complete mass dependence up to order a 2 s 0. Thus in this 
region we will replace data by the theoretical prediction for R{s) as discussed in Section [|. 

In Tab. we present the results for the moments separated according to the three 
different contributions discussed above. The error of the resonance contribution is due to 
the uncertainties of the input parameters. In the case of the charm threshold contribution 
the uncertainty is dominated by the correlated normalization error of approximately 4.3% 
of the BES data. In the continuum region we varied the input parameters as as stated 
in Eq. ([5]) and the renormalization scale as fi — (3 ± 1) GeV. The errors of the three 
contributions are added quadratically. It is illustrating to compare the composition of the 
experimental error for the different moments. Generally speaking, it is dominated by the 
resonance contribution, specifically by the 7% and 9% uncertainty in the leptonic widths 
of the J/^f and respectively. For the moment with n = 1 and to some extent the one 
with n = 2 the improvement in the cross section measurement due to BES (from about 
10 - 20 % systematic error down to 4.3%) was important. The parametric uncertainties 
(from a s and M c ) and the residual /i-dependence which affect Ai c ° nt are small. The higher 
moments (in fact already for n above two) are increasingly dominated by the resonance 
contributions with their 7% uncertainty. 

We use the results of Tab. S together with Eq. fll~9|) in order to obtain in a first step 



m c (3 GeV). Subsequently the result is transformed to the scale-invariant mass m c {m c ) JTJ] 
including the three-loop coefficients of the renormalization group functions. Both results 
can be found in Tab. f|. The starting value for a s (3 GeV) = 0.254±o;[n4 needed for this step 
is obtained from a s (Mz) = 0.118 ± 0.003 by using the renormalization group equations 
and the matching conditions with four loop accuracy. 

The errors listed in Tab. (| receive contributions from the uncertainties in the exper- 
imental moments and the variation of a s (cf. Eq. (|5|)) and /j, — (3 ± 1) GeV in the 
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ft 


KA cxp,res 


A/fexp,cc 


M cont 

Jvl n 


A/f cxp 






xin( n - 1 ) 


A ±U 


xin^- 1 ) 


1 


0.1114(82) 


0.0313(15) 


0.0638(10) 


0.2065(84) 


2 


0.1096(79) 


0.0174(8) 


0.0142(3) 


0.1412(80) 


3 


0.1094(79) 


0.0099(5) 


0.0042(1) 


0.1234(79) 


4 


0.1105(79) 


0.0057(3) 


0.0014(0) 


0.1175(79) 


5 


0.1126(80) 


0.0033(2) 


0.0005(0) 


0.1164(80) 


6 


0.1155(82) 


0.0020(1) 


0.0002(0) 


0.1176(82) 


7 


0.1190(84) 


0.0012(1) 


0.0001(0) 


0.1202(84) 


8 


0.1230(87) 


0.0007(0) 


0.0000(0) 


0.1237(87) 



Table 3: Experimental moments as defined in Eq. (|i~8|) separated according to the contri- 
butions from the resonances, the charm threshold region and the continuum region above 
y/a = 4.8 GeV. 



n 


1 


2 


3 


4 


m c (3 GeV) 


1.027(30) 


0.994(37) 


0.961(59) 


0.997(67) 


m c (m c ) 


1.304(27) 


1.274(34) 


1.244(54) 


1.277(62) 


n 


5 


6 


7 


8 


m c (3 GeV) 


1.094(110) 


1.184(161) 


1.253(182) 


1.307(191) 


m c (m c ) 


1.366(100) 


1.447(146) 


1.510(165) 


1.558(172) 



Table 4: Results for m c (3 GeV) and m c (m c ) in GeV obtained from Eq. ([T9|). 



coefficients C n which are all added quadratically. It is interesting to note that the un- 
certainty in m c (3 GeV) induced by the experimental moments decreases from 0.028 for 
n = 1 to 0.010 for n = 4 whereas the theoretical uncertainty from the renormalization 
scale increases from 0.001 to 0.064, the one from a s from 0.011 to 0.019. 

The moment with n — 1 is evidently least sensitive to nonperturbative contributions 
from condensates, to the Coulombic higher order effects, the variation of [i and the para- 
metric a s dependence. Hence we adopt 

m c (m c ) = 1.304(27) GeV. (22) 



as our final result. Using the two- and three- loop relation |26], |27|, [28[] between the pole- 
and the MS-mass this corresponds to 

Mc (2-ioop) = 1.514(34) GeV, 

Mc (3-ioop) = i 691(35) GeV. (23) 

Our result agrees within the uncertainties with the recent determinations of m c (m c ) in [|K], 
1T| but is comparatively more precise. 
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Figure 3: m c (m c ) for n = 1,2,3 and 4. For each value of n the results from left to 
right correspond the inclusion of terms of order a] and a 2 s in the coefficients C n (cf. 
Eq. (0)). Note, that for n = 3 and n = 4 the errors can not be determined with the help 
of Eq. ( |19D in those cases where only the two-loop corrections of order a s are included 
into the coefficients C n as the equation cannot be solved for m c (3 GeV). 



In Fig. |] we compare the results for m c (m c ) (and its /i dependence) based on the 
theory moments {n = 1,...,4) evaluated up to 0{a\) and 0(a 2 ), respectively. 

The improved stability with increasing order in a s is evident, and the preference for the 
first moment is clearly visible. 

As an alternative we also present the mass values derived from the ratio of moments 
(Eq. (12(f) ) in Tab. [|. The error is obtained from the quadratic combination of the un- 
certainties induced by the resonances, the integration in the charm threshold region, the 
continuum contribution and the independent variation of a s and /i in the coefficients C n . 
It is significantly larger than the one based on the analysis of the moments. The numbers 
are therefore compatible with Eq. (|22"D but do not improve upon this result. 



4 The bottom quark mass 

The same approach is also applicable to the determination of m^. The coefficients C n are 
listed in Tab. |6|. They determine the theoretical moments through Eq. ([16]) where Q c has 
to be replaced by Qb = —1/3. The experimental results for the moments are listed in 
Tab. 0. The contribution from the resonances include Y(IS') up to T(6S) and are given 
by .M^ xp,res . The treatment of the region between y'i = 11.075 GeV and ^/s = 11.2 GeV 
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n/ (n + 1) 


1/2 


2/3 


3/4 


4/5 


m c (3 GeV) 
m c (m c ) 


0.954(64) 
1.238(59) 


0.876(340) 
1.166(318) 


1.052(153) 
1.327(140) 


1.259(261) 
1.515(236) 


n/ (n + 1) 


5/6 


6/7 


7/8 




m c (3 GeV) 
m c (m c ) 


1.434(339) 
1.671(304) 


1.574(386) 
1.795(344) 


1.679(408) 
1.886(361) 





Table 5: Results for m c (3 GeV) and m c (m c ) in GeV obtained from Eq. (|2 



n 


1 


2 


3 


4 


5 


6 


7 


8 




1.0667 


0.4571 


0.2709 


0.1847 


0.1364 


0.1061 


0.0856 


0.0709 




2.5547 


1.1096 


0.5194 


0.2031 


0.0106 


-0.1158 


-0.2033 


-0.2660 




2.1333 


1.8286 


1.6254 


1.4776 


1.3640 


1.2730 


1.1982 


1.1351 


MX) 


3.1590 


3.2319 


2.0677 


1.2204 


0.7023 


0.4304 


0.3359 


0.3701 


r*(21) 


3.4425 


5.0798 


4.5815 


3.4726 


2.1508 


0.7592 


-0.6426 


-2.0281 




0.0889 


1.9048 


3.3185 


4.4945 


5.5127 


6.4182 


7.2388 


7.9929 



Table 6: Coefficients of the photon polarization function in the MS scheme as defined 
in Eqs. flT3|) and (II). n/ = 5 has been adopted which is appropriate for the bottom 
threshold. 



follows | 25| . We assume a linear raise from zero to the pQCD prediction i2&((11.2 GeV) 2 ) 
and the take the contribution itself as an estimate for the error. Note that the contribution 
to the moments, .M® xp ' lm is negligible small. Above 11.2 GeV we use the prediction from 
pQCD for Rb{s) which results in the moments Al^ ont . 

The uncertainties are completely analogous to the charm quark case. The only differ- 
ence concerns the renormalization scale for which we adopt /x = (10 ± 5) GeV. 

The results for m&(10 GeV) and mb(mb) are listed in Tab. A remarkable consistency 
and stability is observed. For n — 1 the error is dominated by the experimental input. 
For n = 3 we obtain ±0.036 from the experimental input, ±0.025 from a s and ±0.020 
from the variation of \i. 

The sensitivity to the inclusion of higher orders is displayed in Fig. [|. Again a signif- 
icant improvement of the stability of our prediction is observed. As our final result we 
adopt 

m b (m b ) = 4.191(51) GeV, (24) 

well consistent with evaluations based on the analysis of bottonium and high spectral 
moments |29], |3(J, |3l], |32], ^3| (see also |22|])[|]. The MS result of Eq. (^4j) corresponds to a 



3 Our result also agrees with the one of Ref. |54j. However, in |34| the bottom quark mass has been 
determined from relatively large moments (7 < n < 15), a highly disputable treatment of the threshold 
has been employed and no stability for small n has been observed. 
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11 


A/f cxp,res 


A^exp.lin 


Accent 


M cxp 

J L n 




xl Q(2n+l) 


xl g(2n+l) 


xl g(2n+l) 


xl g(2n+l) 


1 


1.508(114) 


0.035(35) 


2.913(21) 


4.456(121) 


2 


1.546(109) 


0.028(28) 


1.182(12) 


2.756(113) 


3 


1.600(106) 


0.022(22) 


0.634(8) 


2.256(108) 


4 


1.671(104) 


0.018(18) 


0.381(5) 


2.070(105) 


5 


1.761(103) 


0.014(14) 


0.244(3) 


2.019(104) 


6 


1.870(104) 


0.012(12) 


0.162(2) 


2.044(104) 


7 


1.998(105) 


0.009(9) 


0.111(2) 


2.119(106) 


8 


2.149(108) 


0.007(7) 


0.078(1) 


2.234(109) 



Table 7: Moments for the bottom quark system. 



4.9 



4.8 



4.7 



> 4.6 - 
u 

O 

X 4.5 



4.4 



4.3 



4.2 



4.1 



Figure 4: m&(m.i,) for n = 1,2,3 and 4. For each value of n the results from left to 
right correspond the inclusion of terms of order a°, a] and a 2 s in the coefficients C n (cf. 



Eq. (13)). Note, that the errors for n = 3 and both the central value and the errors for 
n = 4 can not be determined in those cases where only the two-loop corrections of order 
a s are included into the coefficients C n as the corresponding equation cannot be solved 
for m 6 (10 GeV). 



pole mass of [26, 27, 28 



M, 



(2-loop) 
(3— loop) 



4.651(57) GeV. 
4.819(57) GeV, 



(25) 
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n 


1 


2 


3 


4 


m 6 (10 GeV) 


3.665(60) 


3.651(52) 


3.641(48) 


3.655(77) 


m b {m b ) 


4.205(58) 


4.191(51) 


4.181(47) 


4.195(75) 


n 


5 


6 


7 


8 


m b (W GeV) 


3.720(195) 


3.833(293) 


3.965(347) 


4.089(436) 


m b {m b ) 


4.258(189) 


4.367(285) 


4.494(335) 


4.614(420) 



Table 8: Results for m b (10 GeV) and m b {m b ) in GeV obtained from Eq. (|l9|) as described 
in the text. 



nj (n + 1) 


1/2 


2/3 


3/4 


4/5 


m 6 (10 GeV) 
m b (m b ) 


3.636(52) 
4.177(51) 


3.620(51) 
4.161(50) 


3.694(240) 
4.233(233) 


3.915(523) 
4.446(507) 


n/ (n + 1) 


5/6 


6/7 


7/8 




m 6 (10 GeV) 
m b (m b ) 


4.218(717) 
4.737(692) 


4.523(875) 
5.028(839) 


4.786(1.113) 
5.277(1.063) 





Table 9: Results for m b (10 GeV) and m b (m b ) in GeV obtained from Eq. (|20|) as described 
in the text. 



using 0(a 2 s ) and O(a^) accuracy, respectively. 

5 Conclusions 

Recent experimental data for the total cross section a(e + e~ — > hadrons) have been com- 
pared with the up-to-date theoretical prediction of perturbative QCD for those energies 
where perturbation theory is reliable. The excellent agreement justifies the determination 
of the strong coupling a s from the measurements in the energy region between 2 and 
3.73 GeV below the charm threshold and the region between 4.8 and 10.52 GeV above 
charm and below the bottom threshold. Our result a s (5 GeV) = 0.2351q;o47 corresponds 
to a s {M z ) = 0.124t8;8i4 and serves as a useful cross check in the intermediate energy 
region but is less precise than those from r- or Z-boson decays. 

The direct determination of the short distance MS charm quark mass is performed 
with the help of the precise data from the charm threshold region and the three loop 
evaluation of moments in pQCD. Using low moments the approach is insensitive to the 
Coulombic behaviour of the cross section close to threshold and to nonperturbative con- 
densates. The results based on different moments are quite consistent and the moment 
with n = 1 exhibits the least sensitivity towards the parametric dependence on a s and 
the renormalization scale. As our final result we obtain m c (m c ) = 1.304(27) GeV. The 
same approach when applied to the bottom quark gives m b {mh) = 4.191(51) GeV. These 
values are compatable with but more precise than other recent analyses. 
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